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Beamforming with Imperfect CSI: 

An ADMM Approach^ 

Chao Shen*, Tsung-Hui Chang^ Kun-Yu Wang-^-, Zhengding Qiu*, and Chong-Yung Chi-t 



Abstract 



^ I Multi-cell coordinated beamforming (MCBF), where multiple base stations (BSs) collaborate with 

each other in the beamforming design for mitigating the inter-cell interference, has been a subject drawing 
•/^ , great attention recently. Most MCBF designs assume perfect channel state information (CSI) of mobile 

stations (MSs); however CSI errors are inevitable at the BSs in practice. Assuming elliptically bounded 
CSI errors, this paper studies the robust MCBF design problem that minimizes the weighted sum power 
of BSs subject to worst-case signal-to-interference-plus-noise ratio (SINR) constraints on the MSs. Our 
fj ' goal is to devise a distributed optimization method that can obtain the worst-case robust beamforming 

solutions in a decentralized fashion, with only local CSI used at each BS and little backhaul signaling 
for message exchange between BSs. However, the considered problem is difficult to handle even in the 
^ ' centralized form. We first propose an efficient approximation method in the centralized form, based on the 

semidefinite relaxation (SDR) technique. To obtain the robust beamforming solution in a decentralized 
^— N I fashion, we further propose a distributed robust MCBF algorithm, using a distributed convex optimization 

^S) ' technique known as alternating direction method of multipliers (ADMM). We analytically show the 

r-^ , convergence of the proposed distributed robust MCBF algorithm to the optimal centralized solution 

^D I and its better bandwidth efficiency in backhaul signaUng over the existing dual decomposition based 

algorithms. Simulation results are presented to examine the effectiveness of the proposed SDR method 
and the distributed robust MCBF algorithm. 

j^ ■ Index terms— coordinated multi-point (CoMP), multi-cell processing, beamforming, robust beamform- 
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I. Introduction 

Recently, multi-cell processing, or known as coordinated multi-point (CoMP), has drawn great attention 
because it can provide significant system throughput gains compared to the conventional single-cell 
designs [2], [3]. We consider the scenario where the base stations (BSs) are equipped with multiple 
antennas and the mobile stations (MSs) are equipped with single antenna. Each of the BSs employs 
transmit beamforming to communicate with the MSs within its cell. The BSs in different cells, according 
to the principle of interference coordination [2], collaborate with each other to jointly design the beam 
patterns in order to effectively mitigate the inter-cell interference (ICI). To this end, various multi-cell 
coordinated beamforming (MCBF) designs have been proposed [4]-[6]. Most of the MCBF designs 
assume that the BSs are connected with a control center which knows all the MSs' channel state 
information (CSI) and computes the beamforming solution in a centralized manner. In practical multi-cell 
systems, however, obtaining the MCBF solutions in a decentralized fashion using only local CSI at each 
BS is of central importance, thereby drawing the developments of distributed beamforming design methods 
[4]-[10]. The reasons are that 1) the future wireless systems prefer a flat IP architecture where all BSs 
are directly connected with the core network [11]; 2) if the control center is still employed, a distributed 
optimization method can be used to decouple the original problem into multiple parallel subproblems 
with smaller problem size, thus reducing the required computation power of the control center [12]. By 
exploiting the uplink-downlink duality [13], a distributed optimization method was proposed in [4] for a 
power-minimization based MCBF design problem. Game theory based distributed optimization methods 
were proposed recently in [7] for handling the same problem. In [8], the dual decomposition technique 
[14] was used for developing a distributed optimization method for the problem in [4]. In [10], the idea of 
uplink-downlink duality was used for distributed optimization of a different max-min-fair MCBF design 
problem. 

The efficacy of beamforming designs relies on the assumption that the BSs have the perfect CSI of 
MSs. In practical scenarios, however, the BSs can never have perfect CSI, due to, e.g., imperfect channel 
estimation and finite rate feedback [15]. In the multi-cell scenario, it is even more difficult for the BSs to 
obain reliable inter-cell CSI (i.e., the CSI of MSs that belong to the neighboring cells). In the presence of 
CSI errors, the beamforming designs that assume perfect CSI will suffer from performance degradation, 
and the MSs' quality-of-service (QoS) requirements can no longer be guaranteed. In view of this, robust 
MCBF designs, which take into account the CSI errors, are of great importance. Assuming bounded 
CSI errors, e.g., quantization errors, a robust MCBF design was presented in [16] which optimizes the 
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worst-case symbol minimum mean squared error (MMSE) performance of the MSs. 

The focus of this paper is on the worst-case signal-to-interference -plus-noise ratio (SINR) constrained 
MCBF design problem [1], [17], where the weighted sum power of BSs is minimized subject to constraints 
that guarantee worst-case SINR requirements for the MSs. Our goal is to develop a distributed beam- 
forming optimization algorithm for the worst-case robust formulation; however, the considered problem 
itself is difficult to handle even in the centralized form, due to the fact that each of the worst-case 
SINR constraints in essence corresponds to infinitely many nonconvex constraints. The worst-case robust 
design formulation has been studied in the context of single-cell robust beamforming; see [18], [19]. 
However, the robust MCBF formulation is much more challenging since the associated worst-case SINR 
constraints involve CSI errors not only in the desired signal and intra-cell interference terms, but also in 
the inter-cell interference term. To handle this problem, a convex restrictive approximation formulation 
is proposed in [17] which can provide feasible approximate solutions to the robust MCBF problem. 
Distributed optimization algorithms based on dual decomposition and alternating optimization are also 
presented in [17]. However, due to the reduced feasible set, the approximation formulation in [17] is less 
power efficient than the original problem. 

In this paper, we propose a new convex approximation method for the worst-case SINR constrained 
robust MCBF design problem. Our approach is based on a convex approximation technique known as 
semidefinite relaxation (SDR) [20]. SDR has been used in various transmit beamforming designs; see, 
e.g., [21]. By SDR, we obtain a convex approximation formulation for the worst-case robust design 
problem, which, nevertheless, still involves complicated worst-case constraints. We decompose each of 
the worst-case constraints into several simpler worst-case constraints that can be conveniently handled 
by the S-lemma [22]. In particular, by the S-lemma, each decomposed worst-case constraint can be 
reformulated as a linear matrix inequality (LMI) [19]. The resultant approximation formulation, which 
is a convex semidefinite program (SDP), thus can be efficiently solved by interior-point methods [22]. 
We also identify several conditions under which the proposed SDR method can yield the global optimal 
solution of the original worst-case robust MCBF problem. 

We further develop a distributed optimization algorithm that can solve the proposed SDR approxi- 
mation formulation in a decentralized fashion. While the dual decomposition method used in [8], [17] 
is conceptually applicable, we found that the resultant decomposed problems are numerically unstable 
due to lack of strict convexity. To overcome this problem, we instead consider the so called alternating 
direction method of multipliers (ADMM) [23], [24]. ADMM is an advanced dual decomposition method 
that combines the idea of dual decomposition and the augmented Lagrangian method [25], where the 
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latter is often used for bringing numerical robustness to the dual accent method [24] by adding penalty 
terms for strict convexity of the problem. Therefore, in contrast to the conventional dual decomposition 
method [14], ADMM is more numerically stable and faster in convergence [23]. Based on the principle of 
ADMM, we propose a distributed robust MCBF algorithm that is provably able to converge to the global 
optimum of the centralized problem. In particular, by introducing some slack variables that represent the 
worst-case ICI powers, we decompose the SDR problem in a way such that the amount of messages 
required to be exchanged between BSs is much smaller than by the existing algorithms in [8], [17], thus 
reducing the bandwidth overhead of backhaul signaling. 

The rest of this paper is organized as follows. Section II presents the multi-cell system signal model 
and illustrates the impact of imperfect CSI on the performance of MCBF. The considered worst-case 
SINR constrained robust MCBF design formulation is also presented in that section. In Section III, the 
proposed SDR approximation method and its optimality conditions are presented. Using ADMM, the 
proposed distributed robust MCBF algorithm is presented in Section IV. Section V extends the proposed 
method to a fully coordinated scenario where some cell-edge MSs are served simultaneously by multiple 
BSs. Simulation results that demonstrate the effectiveness of the proposed SDR method and distributed 
robust MCBF algorithm are presented in Section VI. Finally, conclusions are drawn in Section VII 

Notations: C", M" and H" stand for the sets of ?i-dimensional complex and real vectors and complex 
Hermitian matrices, respectively. M" denotes the set of ?i-dimensional nonnegative orthant. Column 
vectors and matrices are written in boldfaced lowercase and uppercase letters, e.g., a and A. I„ denotes 
the n X 71 identity matrix, and denotes an all-zero vector (matrix) with appropriate dimension. The su- 
perscripts (•)^, (•)^ and (•)t represent the transpose, (Hermitian) conjugate transpose and pseudo inverse 
operations, respectively. Rank(A) and Tr(A) represent the rank and trace of matrix A, respectively. 
A ^ (^ 0) means that matrix A is positive semidefinite (positive definite). For vector a, ||a|| denotes 
the Euclidean norm. E{-} denotes the statistical expectation. For a variable anmk, where n G {!,..., N}, 
m G {1, . . . , M} and A; € {1, ... , K}, we denote {anmk}k as the set containing Onmi, • • • , CLnmK', while 
we denote {a„mfc} as the set containing all possible a„mfc^ i-C-. oiii, • • • , oiiKj ai2i, ■ ■ ■ , cinmk- 



II. Signal model and Problem Statement 

This section presents the multi-cell downlink system model and the worst-case robust MCBF design 
problem. 
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A. System Signal Model 

Consider a multi-cell downlink system that consists of Nc cells. Each cell is composed of one BS, 
which is equipped with Nt antennas, and K single-antenna MSs. The N^ BSs are assumed to operate over 
a common frequency band and communicate with their K respective MSs using transmit beamforming. 
The scenario under consideration is that each MS is served by only one BS; extension to the scenario 
where one MS is served by multiple BSs will be discussed in Section V. 

We denote BS„ as the nth BS, and MS„fc as the A;th MS in the nth cell, for all n G A/'c = {1, 2, . . . , Nc} 
and /c G /C = {1, 2, . . . . , K}. Let Snk{t) G C be the information data stream for MS„fc, and Wnk G C^' 
be the associated beamforming vector. The transmit signal by BSn is given by 

K 

Xn.(i) = ^W„fcS„fc(*), (1) 

fc=l 

for n = 1, . . . , Nc- The received signal of MS^a; can be expressed as 



Vnkit) = Y^ h^nk [Y^'^miSmiit) + Znk{t) 
m=l \i=l J 



K N, K 

=^nnk^nkSnk{t) + Y^ h^„fcW„iS™(t) + ^1 XI ^mnk^miSmi{t) + Znk{t), (2) 

ij^k mj^n i=l 

where hmnk G C^' denotes the channel vector from BS^ to MS„fc, and Znk{t) G C is the additive noise 
of MS„fc, which is assumed to have zero mean and variance cj^^ > 0. The term Znk{t) may capture 
both the receiver noise and the interference from the other non-coordinated BSs. In (2), the first term is 
the signal of interest for MS„fc, and the second and third terms are the intra-cell interference and ICI, 
respectively. Assume that s„fc(t) are statistically independent, with zero mean and E{|s„fc(t)p} = 1 for 
all n G Mc and A; G /C, and assume that each MS employs single-user detection. By (2), the SLNR of 
MSnfc is given by 

l-i J/ |2 
OTl\TT3 ft ^l^c tu ^Na \ \"-nnk^nk\ 
SINR^fc (^{W„i, . . .,WmK}m=i^{'^mnk}m=l) = "^ ^7"^ • 0) 

Conventional single-cell beamforming designs [26] are developed mainly for handling the intra-cell 
interference only. To take into consideration the ICI, the following multi-cell coordinated beamforming 
(MCBF) design has been considered [27] 

min yZ^niy] l|w„fcf (4a) 

n=l,...,N, "=1 \k=l / 

s.t. SINR^fc ({w^i, . . . , WmK}Z=i^ {hmnk}^=i) > Ink ^k e JC, n G Mc, (4b) 
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where a^ > is the power priority weight for BS„. As seen, the MCBF design jointly optimizes the 
beamforming vectors of all BSs such that the weighted sum power is minimized while the MSs' SINR 
requirements 7^^ > can be fulfilled. It has been shown that problem (4) can be reformulated as a 
convex second-order cone program (SOCP) [27], which can be efficiently solved via standard convex 
solvers, e.g., SeDuMi [28]. 

B. Worst-Case Robust MCBF Design 

The MCBF design in (4) assumes that the BSs have perfect knowledge of the CSIs. As discussed in 
the introduction, in practice, the BSs inevitably suffer from CSI errors. Let h.mnk S C^', n,m ^ Nc, 
A; G /C be the preassumed CSI at the BSs. Then the true CSI can be expressed as 

hmnfc = hmnfc + f^mnk 'i m,n £ J\fc, k £ IC, (5) 

where emnk £ C^' denotes the CSI error vector associated with the true channel hmnk- Our interest 
lies in the bounded CSI errors. Specifically, we assume that each CSI error vector e^n/c satisfies the 
following elliptic model: 

^mnk^'mnk^mnk — ^i v^) 

where Qmnfc G EI^*, Qlrrmk >~ specifies the size and shape of the ellipsoid. When Qm„fc = (l/^mnfc)-'-^! 
where e^^nfc > 0' (6) reduces to the popular spherical error model ||em„fc|p < e^^.^;, [18]. In the presence 
of CSI errors, the non-robust design (4) cannot guarantee the SINR requirement of MSs and consequently 
outage may occur. The following simulation example motivates the need of robust designs: 

Example: We consider a two-cell system {N^ = 2), with two MSs in each cell {K = 2). Each of the 
BSs is equipped with four antennas {Nt = 4). A set of preassumed CSI {hmnA,} is randomly generated 
following the independent and identically distributed (i.i.d.) complex Gaussian distribution with zero 
mean and unit variance (see Section VI-A for the detailed channel model used in the simulation). Using 
the preassumed CSI, an optimal beamforming solution is obtained by solving the MCBF problem (4), 
with a 20 dB target SINR for all the four MSs (i.e., -)nk = 20 dB for n, A; = 1, 2). To examine the impact 
of the CSI errors, we randomly generate 10,000 sets of CSI errors satisfying the following spherical error 
model 

||e„„A.-|P < 0.01, ||e„„fe|p < 0.04 M mi^n, 

for n = 1, . . . , A'^c and A; = 1, . . . , K, and evaluate the achievable SINR values in (3). Figure 1 displays 
the probability distribution of the achievable SINR values of MSn. As seen, it is with very high chance 
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Fig. 1: Distribution of the achievable SINR values of MSn, when the non-robust MCBF design (4) and 
the robust MCBF design (7) are respectively used. The target SINR value is 20 dB. 

that the achieved SINR is smaller than the target value 20 dB, due to the presence of CSI errors. In the 
worst case, the actual SINR value can be even less than 5 dB. ■ 

Our goal is to design the beamforming vectors such that the SINR requirement 7„fc can be guaranteed 
for all possible CSI errors. This idea can be realized by considering the following worst-case robust 
MCBF design [1], [17]: 



AT, 



K 



mm 

W„fc,fc = l,...,E' 

n=l,...,N, "=1 



E"mE 



|W„fc| 



(7a) 



vfc=l 



S.t. SINR„fc (^{w„i, . . . , VfmK}^=lA^mnk + emnfc}m=l ) > Ink 



(7b) 



In comparison with the non-robust design in (4), the above worst-case robust MCBF design can provide 
guaranteed QoS for the MSs, as illustrated in Fig. 1. Solving the robust design problem (7), however, 
is a challenging task. The reasons are that, firstly, each of the SINR constraints is nonconvex, and 
secondly, there are infinitely many such nonconvex SINR constraints due to the worst-case design 
criterion. The restrictive approximation method presented in [17] is able to yield feasible approximate 
solutions, but is less power efficient due to the reduced problem feasible set. In the next section, we 
propose a new approximation method based on relaxation, by applying the convex optimization based 
semidefinite relaxation (SDR) technique [20]. We will further present conditions under which SDR is 
optimal. Simulation results to be presented in Section VI will show that the proposed SDR method is 
more power efficient than the method in [17]. 
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III. Proposed SDR based Method 

A. Solving (7) by SDR and S-Lemma 

Considering that each of the SINR constraint is nonconvex, we first apply SDR to 'linearize' the 
robust MCBF problem (7). To illustrate this, let us express the objective function of problem (7) as 
X]ri=i '^"Sfe=i ^(^nfc^Sc)' ^'^'^ express each of the worst-case SINR constraints in (7b) as 



( I ^ \ . 

\^"^- i^k I 



N^ / K \ 

^ /_^ y^mnk + ^mnk) I / ,'^mi'^mi I y^mnk + ^mnkj + <7„fc V Gj^j^f^Qmnk^mnk < I, m £ Nc- 

The idea of SDR is to replace each rank-one matrix ^nk^^nk ^y ^ general-rank positive semidefinite 
matrix Wnk h 0, by removing the rank-one constraint on W„fc [20]. After applying SDR to (7), we 
obtain the following problem 

mill Va„ y Tr(W„fc) (8a) 



ni^n \'^— 1 / 

^ f^ZnkQmnkemnk < 1, W G jVc, U G A/'c, /c G /C. (8b) 

Note that the SDR problem (8) is convex, since both the objective function and constraints are linear in 
W„fc. However, the SDR problem (8) is still computationally intractable because it involves an infinite 
number of constraints. Fortunately, the infinitely many constraints can be recast as a finite number of 
convex constraints. 

To show this, we first observe that the left-hand side and right-hand side of the first inequality in (8b) 
involve independent CSI errors. Hence, the worst-case constraint in (8b) for MS„fc can be alternatively 
expressed as 



mm 



y] i „ max (hmnfc + emnfc) y^Wmi (hmnfc + e^^fc) )+<T 



(9) 



m^n ^ ™" ~ \i=l 
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By introducing the slack variable 

l^mnk — max I llmnk~T~^mnk I I / , "mi I "■mnk~T~^mnk I v^UJ 

eS„,Q„„.e„„.<l V J \^^ J V J 

as the worst-case ICI power from BSm to MS„fc, for all m G J\fc\{n}, (9) can be written as 

^1) 



(K \ iVe 
W„fc-VW„i [Knk+ennk) > V] *mnfc + 0-^fc- (^ 

By (10) and (11), the worst-case SINR constraint for MS^^ in (8b) can be decoupled into the following 
Nc worst-case constraints: 

(^ K \ N^ 
W„fc— 2^ W„j I \h.nnk+^nnk] > 2_^ tmnk + 0"„,fc V e^nk^nnk^nnk < 1, (12) 
Ink .,, I ^ ' , 

[h^rik + ^mnk) '^^ mi {^mnk + ^m.nkj < trunk V e^nk^mnkfimnk < 1, m G Mc\{n}. (13) 

The ingredient of reformulating (12) and (13) into finite convex constraints is the S-lemma: 

Lemma 1 (S-lemma [22]) Let (l)i{x.) = x^AjX + bf x + x-^bj + a, for i = 0,1, where Aj E M^\ 
bj G C^* and Ci G M. Suppose that there exits an x G C^' such that (pi (x) < 0. Then the two conditions 
are equivalent: 

(1) (poi'^) ^ for all x satisfying (^i(x) < 0; 

(2) There exists A > such that 



Ao 


bo 


+ A 


Ai 


bi 


bo^ 


CO 




bf 


Cl 



>- 0. 



By applying the above S-lemma, we can equivalently recast (12) and (13) as 

^nk (^i »» nil j=l' I'-mnfc J mi -^nnk) 



^nk ~T~ '^nnk^nnk 



Unfch< 



nnk 



"■nnk^nk "-nnk^ nk"-nnk ^nnk 2^ t^nk <^nk 



m^n 



hO, 



(14) 



^ mnk 1^1 " mi/i=l) f-mnfej ^mnfcj 



K \ 

/ J »* mi I ~r -^rnnk^mnk 
i=l / 

K 



'■rank \ Z^ * * m.i 
\i=l 



K 



K 

/ y »» nii ) "-mnk 
i=l 



n™„i. I 2—i ^ mi I ninfc I Z-/ ^ mi I "-mnk ~r ^-mnk ^mnk 



i=l 



^0, m G 7Vc\{n}, 
(15) 



July 12, 2011 



DRAFT 



10 

where U„fc = ^- W„fc — YliJ^k ^m, and \mnk > are slack variables. In summary, one can reformulate 
(8) as the following problem 

min J]aJ J^T¥(W„fc) (16a) 

S.t. *„fe ({W„j}fLi, {tmnk}m, Knk) h ^ n £ Mc, k £ IC (16b) 

*mnfc ({Wmi}^i,tm„fc, Amnfc) ^ V m G 7Vc\{n}, U £ Mc, k £ )C, (16c) 

W„fc ^ 0, Xmnk > V m, n G TVc, fc G /C. (16d) 

Problem (16) is a convex semidefinite program (SDP) which can be efficiently solved by off-the-shelf 
convex solvers [28]. 

B. Optimality Conditions 

An important aspect of SDR is whether the relaxed problem can yield a rank-one solution, i.e., whether 
the optimal solution {W*^} satisfies W*^ = w*^„(w*^)^ for some vi^*^, G C^', for all n, k. If this 
is true, then {vi^*^.} is an optimal solution of the original robust MCBF problem (7). It, therefore, is 
important to investigate the conditions under which the SDR problem (16) can yield a rank-one solution. 
Some provable conditions are given in the following proposition: 

Proposition 1 Suppose that the SDR problem (16) is feasible. Consider the following three conditions: 
CI) K = 1, i.e., there is only one MS in each cell; 

C2) Qnnk = ooItVj for all n, k, i.e., Gnnk = 0/or all n, k, and thus perfect intra-cell CSI {h,„„fc}; 
C3) For the spherical error model, i.e., ||emnfc|P < ^mnk f^^ ^^^ "^' "■' ^' ^^^ ^^^ error bounds {Smnk} 
satisfy 



^mnk S ^mnk Cmct £nnk ^ \ TJ \^ ') 

for all m,n,k, where {s^^^} are some CSI error bounds under which problem (16) is feasible, 

with /* > denoting the associated optimal objective value. 
If any one of the above three conditions is satisfied, then the SDR problem (16) tnust yield a rank-one 
solution, that is, the optimal solution, denoted by {W*^}, must satisfy 

^nk = KkiKkf for all n, k. 

The proof of Proposition 1 is presented in Appendix A. Case CI) of Proposition 1 shows that the 
global optimum of the robust MCBF problem (7) can be attained by solving the SDR problem (16) 
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whenever K = 1. We should mention that, when K = 1, problem (7) can be handled alternatively by a 
convex conic reformulation approach presented in [17]. For the general case of K > 1, C2) shows that if 
the BSs have channel uncertainty only for inter-cell CSI, i.e., {hmnfc} where m ^ n, and perfectly know 
the intra-cell CSI, i.e., {h„„fc}, then rank-one solutions are guaranteed. If errors occur in both intra-cell 
and inter-cell CSI, C3) states that rank-one solutions can also be obtained by solving the SDR problem 
(16), provided that the CSI errors are sufficiently small. 

For a general setup, it is not known theoretically whether the SDR problem (16) has a rank-one solution. 
If the obtained solution is not of rank one, then additional solution approximation procedure, such as 
the Gaussian randomization method [21], can be used for obtaining a rank-one approximate solution to 
problem (7). Quite surprisingly, we found in our simulation tests (see Section VI-A for the setting), that 
problem (16) always yields rank-one optimal {W*;,}. Hence for these problem instances, we can simply 
perform rank-one decomposition of W*^, = w*^(w*^)^. Investigating the reasons behind would be an 
interesting future research; see [29], [30] for recent endeavors. 

IV. Distributed Robust MCBF Algorithm using ADMM 

Solving the SDR problem (16) calls for a control center which computes the beamforming solutions 
in a centralized manner using all the CSI of MSs. As discussed in the introduction, it is desirable to 
obtain the beamforming solutions in a decentralized fashion using local CSI at each BS, i.e., BS„ uses 
{h.nmi}m,i Only, forn = 1, . . . , Nc- A simple approach would be applying the dual decomposition method 
[14], similar to the works in [8], [17]. However, as will be explained later, the dual decomposition method 
is not suitable for problem (16) due to the fact that the decomposed problems lack strict convexity and 
may be unbounded below. To fix this, we propose the use of alternating direction method of multipliers 
(ADMM) [23], [24] for distributed optimization of the SDR problem (16). In the first subsection, we 
briefly review ADMM. In the second subsection, we show how an effective distributed robust MCBF 
algorithm can be developed following the principle of ADMM. 

A. Review of ADMM 

To illustrate the idea of ADMM, let us consider the following convex optimization problem [23], [24] 

mill F(x) + G(z) (I8a) 

xeK",z6K'" 

s.t. xe5i, z e52, (18b) 

z = Ax, (I8c) 
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where F : M" i— )• M and G : M™ i— )• M are convex functions, A is an rii x n matrix, and Si C M" and 
^2 C M™ are nonempty polyhedral sets. Assume that problem (18) is solvable and strong duality holds. 
ADMM considers the following penalty augmented problem 

mill F(x) + G(z) + ^||Ax-z|P (19a) 

xGR",zGK'" Z 

s.t. X G 5i, z G 52, (19b) 

z = Ax, (19c) 

where c > is the penalty parameter. It is easy to see that (19) is essentially equivalent to (18) owing 
to (19c). The penalty term |||Ax — z|p brings strict convexity; as seen, problem (19) is strictly convex 
with respect to either x or z. 

The second ingredient of ADMM is dual decomposition [24] where the dual problem of (19) is 
concerned: 



max < 



mill F(x) + G(z) + f||Ax-zf + ^^(Ax 



'■ • (20) 



s.t. X G 5i, z G ^2 

in which ^ G W^ is the dual variable associated with the constraint (19c). Given a dual variable ^, the 
inner problem is a convex problem and can be efficiently solved. The outer variable ^ can be updated by 
the subgradient method [22]. In a standard dual optimization procedure, one usually updates the outer 
variable ^ when the associated inner problem has been solved with the global optimum. For example, 
one can use the nonlinear Gauss-Seidel method [23] to optimally solve the inner problem. Specifically, 
in the nonlinear Gauss-Seidel method, one iteratively solves the following two subproblems 

z(g + 1) = arg min |g(z) - i{q)^z + ^||Ax(g) - zfj , (21a) 

x{q + 1) = arg mill |f(x) + ^{qfAx + ^||Ax - z(g + l)f 1 (21b) 

xe<Si L 2 J 

until convergence, where q is the iteration number. Instead, ADMM, as its name suggests, alternatively 

performs one iteration of the Gauss-Seidel step (21) and one step of the outer subgradient update for 

speeding up its convergence. The steps of ADMM are summarized in Algorithm 1. 

ADMM can actually converge to the global optimum of problem (18) under relatively loose conditions; 

see the following lemma from [23, Proposition 4.2]: 

Lemma 2 [23, Proposition 4.2] Assume that Si is bounded or that A^A is invertible. A sequence 
{x(g),z(g),^((7)} generated by Algorithm 1 is bounded, and every limit point of {x{q),z{q)} is an 
optimal solution of the original problem (18). 
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Algorithm 1 ADMM 



1: Set q = 0, choose c > 0, 

2: Initializes ^{q) and x(g); 

3: repeat 

4: z(g + l) = argmin{G(z) -$,^{q)z + §||Ax(g) -z|p} 

ZGO2 

5: x(g + 1) = argmin{F(x) + ^^(g)Ax + §||Ax - z(g + l)|p} 

xGoi 

6: ^(<?+l)= ^(9)+c(Ax(g + l)-z(g + l)) 

7: g:=g+l; 

8: until the predefined stopping criterion is satisfied. 



B. Applying ADMM to Problem (16) 

Applying ADMM to our SDR problem (16) for distributed optimization is not as straightforward 
because the constraints in (16) are intricately coupled. Our intension in this subsection is to reformulate 
problem (16) such that the corresponding ADMM steps in Algorithm 1 are decomposable and thus 
problem (16) can be solved in a distributed fashion. 

To this end, we first introduce the following two auxiliary variables 

K N, 

Pn = ^Tr(W„fc), Tnk = ^ tmnk (22) 

fc=l m^n 

for all n G J\fc and k ^ K,, where p„ represents the transmission power of BS„, and T^^ stands for the 
total worst-case ICI power from the neighboring BSs to MS„fc. Then problem (16) can be rewritten as 

AT, 

min V" UnPn (23a) 

{W„fchO}, ^ 

{A„„,>0},{t„„fc} "-^ 
{p.},{T„.} 

S.t. *„fc {{^ni}f=l,Tnk, Xnnk) ^ V n G AA^, /c G /C, (23b) 

*mnfc {{^mi}f=l,tmnk, ^mnk) ^0, M m £ Nc\{n}, 71 £ Mc, k £ IC, (23c) 

K 

Y, Tr(W„fc) = p„ V n G AA„ (23d) 

fc=i 

y^ trnnk = T^k ^ n £ Mc, k £ /C. (23e) 



It is interesting to observe from (23b) that each MSnk concerns only the total worst-case ICI power Tnk 
instead of the individual worst-case ICI powers {tmnk}- As will be clear later, introducing such slack 
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variables {Tnk} will reduce the backhaul signaling overhead for distributed optimization. 

Note that, without loss of generality, we can interchange the subindices m and n in (23c). Hence, the 
constraints in (23b) to (23d) can be decomposed into Nc independent convex sets: 



^n \\ \^^ nk Sk: \^nmkSm,ki X-'-nkJki X'^nmk Sm,kTPn 

*n,fc ({Wni}£i,7;,fc,A„„fc) ^ V fc G /C, 

*nmfc {V^ni}fLl,tnmk,Kmk) ^ V m G A/'c \ {n}, k £ JC, 

Kmk > V m G TVc, fc G /C, 

E TliWnk) =Pn, WnkhO^kelC 
k=l 



(24) 



for n = 1, 



, Nc- Further define the following variables 



tl21,- ■ ■ ,h2K, •••, *7V,(Are-l)l' • • • !%c(Afc-l)-ft' 



.Nc{Nc-l)K 



tnl, 



> -^riK ■) '■nil) 



<tnlK-,tn2li 



, tnN^K 



^NcK 



, n = l,...,Nc, 



(25a) 
(25b) 



where t collects all the ICI variables, and t„ collects variables {Tnk}k=i ^^id {tnmk}m,k (where m ^ n) 
that are relevant only to BS„. It is not difficult to check that there exists a linear mapping matrix 
E„ G ]K^<:^^x^<:(^<:-i)^^ wMch contalns only elements either equal to one or zero, such that 

tn = E„t, (26) 

for all n = 1, ... , N^. By (24), (25a) and (26), we can rewrite problem (23) in a compact form as 



JVc 



mm 

{W„.},{A_4,{t„},{p„} 



,t ^ 



n=\ 



S.t. W ^ nksk-i l-^nm/c Jm,fc) ^n^Vn) G ^ni '^ — -Lj • • • j -'^C! 

tri = E„t, n = 1,... ,iVc- 



(27a) 

(27b) 
(27c) 



Before applying ADMM, let us first see why the conventional dual decomposition method [14] is not 
suitable for problem (27). Consider the dual problem of (27) 



max < mm 

n=l,...,Afe V n=l A^c 






(28) 



_ _ n.=l n=l 

where Vn G M^^^, n = 1, . . . ^N^, are the dual variables associated with constraints (27c). While the 



inner minimization problem of (28) is obviously decomposable, given i^i , . . . , iv^v^ , it is possible for 

— oo, i.e., the inner minimization problem is 



one to obtain an inner solution of t„ such that — i/Jt„ 
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unbounded below, due to the unboundedness of the feasible sets C„ (see (24), (14) and (15)). In fact, our 
numerical experience shows that this undesired situation happens very often, especially when Nc > 2. 

To overcome this issue, we apply the augmented Lagrangian method to (27) according to the principle 
of ADMM; this leads to the following problem: 



JV, 



c sr-^ ,,„ ,,9 c 






n=\ 



2 A^ii-.- -.M, . 2 

n=l n=l 



S.t. ({W„fc}fc,{A„,mA:}m,fc,tn,Pn) G C„, n = 1, . . . , A^c, 



t„ = E„t,n= l,...,iVc, 



Vn = Pn-,n = \,...,Nc 



(29a) 

(29b) 
(29c) 
(29d) 



where /9„ > 0, n = 1, . . . , A^c> are slack variables, which are introduced in order to impose the penalty 
term | X]n=i ^Pn — Pn) ■ Problem (29) is equivalent to problem (27), but the added penalty terms can 
resolve the numerically unbound below issue; see [23]. 

Now we are ready to apply ADMM. Consider the following correspondences between (29) and (19): 



X=[t ,pi,...,pivj , Z = [ti,...,t^^,pi,...,pArJ , A 



E 

Itv. 



(30) 



N, 



F(x)^0, G(z)^ J]a„p„, 

n=l 

c A (rNciN^-l)K+N, 

^2 = {[ti , . . . ,t^^,pi, . . . ,PArJ |({W„fc}fc,{Ar!,mfc}m,fc,tn,p„) G C„,n = 1, . 
/. A r T T iT 

i= [l/i ,...,iv^ ,/ii,...,/i7vJ , 



•,iVc}, 



where E = E?^ • • • ET^ , and Un G M^^^, /i„ G M, n = 1, . . . , Nc, are the dual variables associated 
with constraints (29c) and (29d), respectively. According to Algorithm 1, the corresponding ADMM step 
4 for problem (29) is to solve the following problem 



min Y] [anPn + - ||Ent(g) - t„||^ + - {pn{q) - Pnf " i^n l^)*" " P'n{q)Pn 

{W„fc}fc,{A„„fc}„,fc,t„,p. — \ -* '' 

n=l,...,N, 



n=l 



2 11 - v./ .M, 2 

S.t. (\ VV jj/j j-fc, -[Ajjf^fcJ-n^jfc) I'n) Pnj G Cfi, ?Z ^ i, . . . , iVc. 



(31) 
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As one can see, problem (31) can be decomposed as the following Nc problems: 

{tn{q + l),Pn{q + l)} = 

arg min a„p„ + - ||E„t(g) - t„f + - (p„,(g) -p„)^ - i/J(g)t„, - /i„(g)p„, 

S.t. \\^ nksk^X^nmklm^k^^mVn) & '--n) \~'^) 

for n = 1, . . . , A'^c- Since problem (32) is convex, it can be efficiently solved. 

Secondly, the corresponding ADMM step 5 is given by solving the following two problems: 

t(g + l) = arg min ^J]||E„t-t„(g + l)f + J]i.^(g)E„t, (33) 

{pn{q + l)}^=i = arg min - V(pn - Pn{q + 1))^ + V fin{q)Pn- (34) 

p„>0, 2 ^ ^ 

n=l,...,Af, "=1 "=1 

Because both (33) and (34) are convex quadratic problems, they have closed-form solutions given by 

t{q+l)=B^(i{q+l)--^i>{q)Y (35a) 

Pn{q + 1) = Pn{q + 1) fJ-niq), H = 1, . . . , Nc, (35b) 

c 

where i{q + 1) = [tf (g + 1), . . . , t%^{q + 1)]^ and i>{q) = [i^fiq), ..., i^%^iq)f . 

Finally, the corresponding ADMM step 6 is given by the following dual variable update 

i^niq + 1) = Mq) + c (E„t(g + 1) - t„(g + 1)) , (36a) 

IJn{q + 1) = Mn(5) + C (p„(g + 1) - Pn(5 + 1)) , (36b) 

for n = 1, . . . , Nc- 

It is important to note that the ADMM steps (32), (35) and (36) can be implemented in a distributed 
fashion. Essentially, given the knowledge of local CSI {hnmfc}m,fc. the optimization problem (32) can 
be independently solved by BS,i, for all n = 1, . . . , Nc- After that, each BS^ broadcasts its new t„ to 
the other BSs. With the knowledge of {t„}, each BS can compute the public variable t by step (35a). 
Moreover, Pn{q + 1). '^n{q + 1) and Pn{q + 1) in (35b), (36a) and (36b) all can be independently updated 
by each BS.„, for n = 1, . . . ,N(.. Summarizing the above steps, we thus obtain the distributed robust 
MCBF algorithm in Algorithm 2. 

Algorithm 2 is guaranteed to converge to the global optimum of the SDR problem (16). Specifically, 
one can verify that the matrix A in (30) satisfies 

E:=iE^e„ 

1n^ 



A^A 



>- 0. (37) 
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Algorithm 2 Proposed Distributed Robust MCBF Algorithm: 



1: Input a set of initial variables {^'n(0),/in(0), t(0),p„(0)}„JL]^ that are known to all BSs; choose a 

penalty parameter c > 0. 
2: Set q = 0. 
3: repeat 
4: Each BS„ solves the local beamforming design problem (32) to obtain the local ICI variables 

t„((/+l) and the local power pn{q+l). 
5: Each BS„ broadcasts its local ICI variable t„ to the other BSs, e.g., via the backhaul network. 
6: Each BS„ updates the public ICI t and p„ by (35a) and (35b), respectively. 
7: Each BS„ updates the dual variable Un an ^„ by (36a) and (36b), respectively. 
8: Setg:=g'+1; 
9: until the predefined stopping criterion is met. 

By Lemma 2 and by the fact of (37), we obtain the following result on the convergence of Algorithm 2: 

Proposition 2 Consider the proposed distributed robust MCBF algorithm in Algorithm 2. The iter- 
ates t{q), {pn{(l),'^n{Q), Pn{Q)}n=i '^^^ {^n{(l) , fJ-niq)} n=i ^'^^ respectively converge to the optimal 
primal and dual solutions of problem (29) as q ^ oo. When the algorithm converges, the optimal 
{W„i, . . . ,'WnK}n=i obtained in Step 4 is a global optimal solution of the SDR problem (16). 

Three remarks regarding the proposed distributed robust MCBF algorithm are in order. 

Remark 1 In Algorithm 2, each BS„ needs to exchange its local ICI power vector t„(g + 1) with the 

k=l 



Other BSs through the backhaul network. Since t„(g+l) contains K total incoming ICI powers {Tnk} ' 
and {Nc — l)K outgoing ICI powers {tnmk}m^n,k, the total backhaul signaling is N^K real variables for 
each iteration. By contrast, for the distributed algorithms in [8], [17], each BS has to exchange {Nc — l)K 
real variables for incoming ICI powers and {Nc — l)K real variables for outgoing ICI powers, and thus 
a total number of 2Nc{Nc — l)K real variables need to be exchanged for each iteration. For Nc > 2, the 
proposed algorithm is clearly more backhaul bandwidth efficient. For example, if Nc = 6, the required 
backhaul signaling of the proposed algorithm is about 60% of that of the algorithms in [8], [17]. 

Remark 2 Interestingly, Algorithm 2 can be interpreted as an adaptive ICI regularization scheme where 
the cooperative BSs gradually control their own beamforming solutions until a consensus on the induced 
ICI powers among BSs is reached. To further explain it. Step 6 of the algorithm can be regarded as a step 
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that computes the tentatively consentient ICI power vector t(g + 1) based on the locally optimized ICI 
power vectors {tn{q)}n=i- I^i Step 7, the BSs then update the dual variables according to the difference 
between the consentient ICI powers and local ICI powers. Once the algorithm converges (which implies 
E„t((7 + 1) = tn{q + 1) for all n), all the BSs achieve a global consensus on the ICI powers, and hence 
the beamforming solutions obtained in Step 4 are globally optimal (by Proposition 2). 

Remark 3 Since ADMM operates in the dual domain, the obtained {Wnk} and {Xnmk} in Step 4 may 
not be feasible to the primal SDR problem (16). To fix this, each BS may perform one more primal 
optimization of 

K 

Yl Tr(W„fc) (38) 



mm 

{w„fc>:O}fc,{A„™fc>0}„ 



k=l 



S.t. *„fc ({W„i}j^i, {tmnk{q + l)}m5^n, A„„fc) ^ V fc G /C, 

^nmk {{^ ni}f=l:tnmk{q + 1), A„mfc) ^ V ?n G A^c \ {n}, fc G /C, 

using the tentatively consented ICI power vector t{q + 1). The obtained {Wnk} and {Xnmk} then must 
be feasible to the SDR problem (16), provided that problem (38) is feasible for all BSs. If at least one of 
the BSs declares infeasibiUty of (38), then more iterations are needed for Algorithm 2 since the algorithm 
may stop too early to reach a reasonable consensus on the global ICI t{q + 1). 

V. Extension to Fully Coordinated BSs 

In this section, we extend the robust MCBF design to the scenario where some of the MSs are served 
simultaneously by multiple BSs. The scenario may occur, for example, when some of the MSs are near the 
cell boundary and thus desire to receive the information signal sent from multiple BSs for the guaranteed 
QoS. To simultaneously serve these MSs, the BSs have to be fully coordinated, with shared data streams 
and CSI of these cell-edge MSs [2]. Assume that there are L such cell-edges MSs, in addition to the 
NcK intra-cell MSs that are served solely by their respective BSs. The transmit signal of BS„ is given 
by 

L 

i„,(t) = x„(t) + Y, ^nideit), (39) 

e=i 

where x„(t) is defined in (1) which is intended for the K intra-cell MSs, d£{t) is the data stream for 
the ith cell-edge MS, and {^e £ C^' is the beamforming vector of BS„ for sending d£{t). The received 
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signals of intra-cell MS„fc and cell-edge MS i are respectively given by 

ylT\t) = E hL.x.(t - T^r^) + E E hL.fn,rf,(t - rir^) + Zn,{t\ (40) 

»n=l m=l j=l 

m=l m=l j=l 

for n G A^c, A; G /C and £ G £ = {1,...,L}, where g^e G C^* is the channel vector from BSm to 
cell-edge MS i, and ze{t) is the background noise at cell-edge MS £, which is assumed to be zero mean 
and with variance aj > 0. Note from (40) and (41) that we have taken into account the inevitable time 
delays r^^^r^^™^ > between the BSs and MSs [6], [31]. Assume that r^^"^ / r^^"^ for all m / n, 
and that each di{t) is temporally uncorrelated with zero mean and unit variance. The receiver SINRs 
corresponding to (40) and (41) are given by [6] 

(Intra) fj^^^ . rn ^ r, , -.AT, 






(42) 



K ^ N, K ^ N, L 

Y. \^nnk'^rii\ + Yl Y \^mnk'^rn.i\ + Yl Y l^mnJrnjl'^ + Cr^^k 
i=ik m^n 1=1 m=l j=l 

, , Y \Sm£^me\ 

SmRf"^''^ ({w^„}, {f^,}, {g™4^Li) = N. K ^ "^;' L N. 1 -■ (43) 

Y Y iSme'^rnil + Y Y ISm^'^mjl + ^£ 
m=li=l j^£m=l 

Our goal here is, again, to find the beamforming vectors that are robust against the possible CSI errors. 
As the channel error model for intra-cell MSs, we model the cell-edge MSs' channel as 

gmi = Emi + Vm£ V m G A/;, ^ G £, (44) 

where gmi G C^' is the preassumed CSI, and Vmi G C^' is the CSI error satisfying v^^Q^fV^^ < 1 
in which Qmi >- 0. We consider the following worst-case robust formulation: 

N, / K L \ 

.(Intra) ( r 1 rf 1 rr. , ^ \N^ 



s.t. SINR^^. '^ (^{w™},{fmj},{hm„,fc + em„,fc}m'=ij > 7nfc 

V e^nk^rnnkernnk < 1, m e Mc, n e Mc, fc G /C, (45b) 
gj^j^(Edge) J^^^^.^^ ^^^^.^^ ^ .^^ ^ V„4^Li) > 7^ 

V ^mAmi^mi < 1, m G A/;, £ G /:. (45c) 
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The proposed method based on SDR and S-lemma in Section III-A can be used to handle the above 
problem (45) as well. Firstly, replace each w„fcw^^ and each ^ndni ^y general-rank W„fc >z and 
Fn^ ^ 0, respectively. Secondly, follow the steps as in (8a) to (15) to decouple and transform the worst- 
case constraints in (45b) and (45c) into a finite number of LMIs. The obtained SDR problem can be 
shown to be the following SDP: 



N^ / K L 

mill J] an J] Tr(W„fc) + J^ Tv{¥. 

{F„.}, 

{A„a,K„a 



S.t. *Pnk \\^nisi=l7\'^mnksm}'^nnk) 



nij 



^ mnk \\^ Tniii=\^'^mnk) ^ 



mnk] 



W„fc ^ 0, Xmnk >Oym,n€Mc, k £ IC, 



m=l 




Fmi hO, A„f > V ?n G J\fc, i G C, 






H 



+ 



(46) 




nnk 



mnk 



H 



^0, 



H 



y y m ^ n, 



XmeQme 

—^mi—lime 



hO, 



where {Xmi}, {Vmi} are slack variables. For the spherical error model, a sufficient condition for the 
tightness of SDR, which is similar to Proposition 1, can be shown to 



^nmk — ^nmki ^n£ — ^n£ 



£nnk < \ ^^^, and Eni < \ ^ V n, m, k, £, 



(47) 



where {Snmk}, {^ne} are the CSI error bounds for which problem (46) is feasible, and 5* > is the 
associated optimal objective value. The condition in (47) implies that problem (46) can attain the global 
optimum of (45) if the CSI errors are sufficiently small. 

A distributed optimization algorithm for problem (46) can also be developed by applying ADMM, using 
the same ideas as presented in Section IV-B for problem (16). We will provide a simulation example in 
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Section VI to demonstrate the efficacy of formulation (46) in providing guaranteed QoS for the cell-edge 
MSs. 

VI. Simulation results 

In this section, some simulation results are presented to examine the performance of the proposed 
robust MCBF design and the distributed optimization algorithm (Algorithm 2). The performance of the 
robust fully coordinated BF design in the previous section will also be examined. 

A. Simulation Setting 

In the simulations, we not only consider the small scale channel fading but also the large scale fading 
effects such as shadowing and path loss, in order to simulate the multi-cell scenario. Specifically, we 
follow the channel model [4], [32]: 

"■mnk — -L" ' Ymnk ' fmnk ' \"-mnk ~r ^mnkji x^o) 

where the exponential factor is due to the path loss depending on the distance between the mth BS and 
MSnfc (denoted by dmnk in kilometers), if^mnk reflects the shadowing effect, and i^mnk represents the 
transmit-receive antenna gain. The term inside the parentheses in (48) denotes the small scale fading 
which consists of the preassumed CSI h.mnk and the CSI error e^nk- As seen from (48), it is assumed 
that the BSs can accurately track the large scale fading, and suffers only from the small scale CSI errors. 
The inter-BS distance is 500 meters, and the locations of the MSs in each cell are randomly determined 
with the distance to the serving BS at least 35 meters, i.e., dnnk ^ 0.035 for all n, k. The shadowing 
coefficient ijjmnk follows the log-normal distribution with zero mean and standard deviation equal to 8. 
The elements of the preassumed CSI {hmnfc} are i.i.d. complex Gaussian random variables with zero 
mean and unit variance. We also assume that all MSs have the same noise power spectral density equal 
to -162 dBm/Hz (-92 dBm over a 10 MHz bandwidth), and each BS has a maximum power limit 46 
dBm [32]. The SINR requirements of MSs are set the same, i.e., ^mnk — 7> and each link has the same 
antenna gain f^mnk = 15 dBi. The power weight a„ for BS„ is set to one for all n (i.e., sum power). 
For the CSI errors, the spherical error model is considered, i.e., Qmnfc = (l/^mnfc)-'-^* ^'^^ ^^^ "^' ^ ^'^'^ 
k. If not mentioned specifically, the error radii e^nk are set the same and denoted by e. 

B. Performance Comparison with Existing Methods 

For the robust MCBF design (7), we first compare the proposed SDR method with the convex 
restrictive approximation method in [17]. The single-cell beamforming (SCBF) design with independent 
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O Robust MCBF (Proposed SDR method) 
-0— Robust MCBF (Method in [17]) 
-H — Robust SCBF 




O Robust MCBF (Proposed SDR method) 
-0— Robust MCBF (Method in [17]) 
-B— Robust SCBF 



■t 




(b)iVe 



Fig. 2: Feasibility rate (%) versus SINR requirement ^ ioi K = A, Nt = 6, e = 0.1. 

ICI constraints [1], [33] is also compared. All the design formulations are solved by SeDuMi [28]. 

We first present the feasibility rates of the three beamforming designs. We say that the formulation 
under test is feasible if it can yield an optimal solution with each BS's power no greater than 46 dBm. 
Figure 2 presents the simulation results ior K = 4 (MSs/cell), Nt = 6 (antennas/BS), and (a) N^ = 2 and 
(b) Nc = 3, respectively. The CSI error radius e is set to 0.1. Seven thousand channel realizations are 
tested. We can observe from this figure that the robust MCBF design exhibits a much higher feasibility 
rate than the SCBF design, showing the improved capability of coordinated beamforming by exploiting 
the degrees of freedom provided by multiple BSs. Secondly, we can see that the robust MCBF design 
using the proposed SDR method exhibits a slightly higher feasibility rate than that using the method in 
[17]. We should emphasize that in the simulation tests, the SDR problem (8) all yields rank-one solutions. 
Hence, the feasibility rate of the SDR method is in fact that of the original problem (7). 

Next, we examine the average transmission sum powers of the robust MCBF design. As a performance 
benchmark, we also present the average sum power of the non-robust MCBF design (4). Figure 3(a) shows 
the results of average sum power (dBm) versus the SINR requirement 7 for Nc = 2, K = A, Nt = 6, 
and e = 0.1. 7,000 channel realizations are tested and each of the results in Figure 3(a) is obtained by 
averaging over the feasible channel realizations at each 7. One can observe from the figure that, as a 
price for worst-case performance guarantee, the robust MCBF designs require higher average transmission 
powers than the non-robust design. Comparing the proposed SDR method with the method in [17], one 
can see that the proposed SDR method is more power efficient. For example, for 7 = 10 dB, the proposed 
SDR method consumes around 24 dBm while the method in [17] requires 29 dBm. Figure 3(b) displays 
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(a) e = 0.1 (b) 7 = 10 dB 

Fig. 3: Average transmission sum power (dBm) of various methods for Nc = 2, K = 4: and Nt = 6. 



the results of average transmission sum power (dBm) versus the CSI error radius e, for 7 = 10 dB. As 
seen, the proposed SDR method is much more power efficient than the method in [17]. 

C. Performance of Proposed Distributed Robust MCBF Algorithm 

Now, let us examine the performance of the proposed distributed robust MCBF algorithm (Algorithm 
2). In the simulations, the initial input values {^'n(0),/^n(0),t(0),p„(0)}„^]^ are all set to zero. The 
augmented penalty parameter c is not fixed as a constant; instead it follows the rule 

{qc{q) if c(g) < 1, 



c{q + l) 



(49) 



[ 1 otherwise, 

where c(0) = 10^^. As one can verify that c{q) will reach the value of one after nine iterations. 
Algorithm 2 following (49) will still converge to the global optimum for a sufficiently large q, according to 
Proposition 2. Firstly, we compare the optimal sum power of the centralized problem (7) with that obtained 
by Algorithm 2. The simulation results by testing over 50 randomly generated channel realizations are 
presented in Fig. 4, under various simulation settings. From Fig. 4(a) and Fig. 4(b), where Nc = 2 and 
K = 2 and K = 4, respectively, we can observe that Algorithm 2 can yield near-optimal solutions within 
50 iterations. As observed, for most of the cases, 10 and 20 iterations are quite sufficient for the scenarios 
in Fig. 4(a) and Fig. 4(b), respectively. When the the number of cells Nc increases to three {Nc = 3), 
as shown in Fig. 4(c), 25 iterations are sufficient to obtain a near-optimal solution. General speaking, as 
the number of cells and that of MSs increase, the number of iterations needed to achieve a near-optimal 
performance also increases. As seen from Fig. 4(d), where Nc = 8, at leaf 100 iterations are required. 
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(d) Nc = ?,,K = l,Nt=A 
Fig. 4: Sum power comparison between the centralized robust MCBF solution and the distributed robust 
MCBF solution (by Algorithm 2) over 50 randomly generated channel realizations, for 7 = 10 dB and 

e = 0.05. 

To further look into the convergence behavior of Algorithm 2, we show in Fig. 5(a) the typical 
convergence curves of Algorithm 2 in the scenarios considered in Fig. 4(a) to Fig. 4(c). In Fig. 5, 
the normalized power accuracy is defined as 

\P*{q)-P*\ 



Normalized power accuracy 



p^ 



(50) 



^N, 



where P*{q) = X]n=i^«('?) ^^ ^^^ ^^™ power at iteration q, and P* denotes the centralized solution of 
problem (16). We can see from Fig. 5(a) that Algorithm 2 can yield a solution with the normalized power 
accuracy smaller than 0.1 within 50 iterations. Figure 5(b) presents the convergence curves of Algorithm 
2 for Nc = 2, K = 2, Nt = 4, and various CSI error radii and SINR requirements. It can be seen from 
this figure that the convergence speed can be slowed down as the CSI error radius or SINR requirement 
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increases. Nevertheless, for the scenarios considered in Fig. 5(b), less than 40 iterations are needed for 
achieving 0.01 normalized power accuracy. The simulation results presented in Fig. 4 and Fig.5 well 
demonstrate the convergence of Algorithm 2, as stated in Proposition 2. 
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Fig. 5: Typical convergence curves of Algorithm 2 under various simulation settings. 

D. Performance of Robust Fully Coordinated BF 

In this subsection, we examine the effectiveness of the robust fully coordinated BF design (45) in 
serving the cell-edge MSs. To this end, let us consider a three-cell system (Nc = 3) with two MSs in 
each cell (K = 2), as illustrated in Fig. 6. We divide each cell into two parts, namely, the intra-cell region 
and the cell-edge region. In particular, the inter-BS distance is set to 500 meters and the radius for the 
intra-cell region is 235 meters. In each cell, the position of one of the MSs is randomly chosen within the 
intra-cell region; while the other MS is randomly located in the cell-edge region within the equilateral 
triangle formed by the three BSs (see Fig. 6). As the robust fully coordinated BF design (45) is applied, the 
three MSs in the cell-edge regions will be served simultaneously by the three BSs, i.e., K = 1 and L = 3 
in (45). Figure 7 shows the performance comparison results of the robust fully coordinated BF design 

(45) and the robust MCBF design (7) by testing over 17,000 channel reahzations. The SDR formulation 

(46) is used as an approximation to (45). It is found in this simulation test that SDR formulation (46) 
always yields rank-one solutions; hence the obtained solution is exactly the optimal solution of (45) for 
the tested problem instances. From Fig. 7, we can observe that the robust fully coordinated BF design 
is more feasible and is more power efficient (for around 3 dB) than the robust MCBF design in serving 
cell-edge MSs. 
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Fig. 6: Illustration of the simulation scenario for the robust fully coordinated BF design (45). 




(a) (b) 

Fig. 7: Performance comparison results of the robust MCBF design (7) and the fully coordinated BF 
design (45), for Nc = 2,, K = 2, Nt = 4, and e = 0.1. 

VII. Conclusions 

In this paper, we have investigated the worst-case SINR constrained robust MCBF design problem [in 
(7)]. While the robust design problem involves complicated nonconvex worst-case SINR constraints, we 
have presented an efficient approximation method based on SDR. We have shown that when there is 
only one MS in each cell or when the CSI errors are sufficiently small, the proposed SDR method can 
yield the global optimal solution to the original problem (Proposition 1). Moreover, by using ADMM, we 
have proposed a distributed robust MCBF algorithm (Algorithm 2). The proposed distributed algorithm 
is appealing because it is proven to converge to the global optimum of the centralized problem, with 
a much less backhaul signaling overhead compared to the existing methods (Proposition 2). Extension 
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of the proposed SDR method to a fully coordinated scenario has been presented as well. The presented 
simulation results have shown that the proposed SDR method is more power efficient than the existing 
method, and that the proposed distributed optimization algorithm can obtain beamforming solutions with 
the normalized power accuracy smaller than 0.1 in twenties of iterations for the typical scenario of 
Nc < 3. The robust fully coordinated BF design (45) has also been shown more power efficient if there 
are MSs located in the cell boundary. 



Appendix A 
Proof of Proposition 1 

Proof of case CI): We rewrite problem (16) for i^ = 1 as follows 



N, 



mm 

{w 

{A„„},{t,„„},„^ 



in Yl (^nTr{Wn) 

"'\ n=l 



S.t. ^n \^ m X'^ninSm^m '^nn) Cl "; 

^ nin (, '' m j i^nin i '^mn ) CI " " TTl ^ Tl, 

W„ y 0, Xmn > V m, n G Mc, 



(A. la) 

(A. lb) 
(A.lc) 
(A. Id) 



where the subindices k and i are removed for notational simplicity. Proposition 1 can be proved by 
investigating the KKT conditions of problem (A. 1). According to the KKT conditions, one can verify 
that the optimal {W*}, {AJ^„},{t^„} satisfy the following conditions 



Z*W* = W* ^ 



IJV, + 



E [^N, h 
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^nn > 0) 
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(A.4) 
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'^nn '-l^n [^ n^ X^mnJrn^m \ 



(A.9) 



and Y*^ h 0, m ^ n, are the associated dual variables. Specifically, to show (A.7), we note that if 
Xnn = 0, then it follows from (14) that 

1 

which contradicts with (A. lb). 

The first step of the proof is to show that W* has rank one whenever Y*„ has rank one. Suppose 
that Y* 



m^n 



E 






yy^ where y G C^*+\ i.e., Y*„ is of rank one. Let X„ = X.II ^J = In, + 



^Nt "-nm 



\H 



y 0. Then 
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Rank(Z*) = Rank | X^^x^/^ _ _ 
Rank I In, - X'^/^ 
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1 



7n. 



^Nt ^nn 
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yy^ 


"■nn 


yy^ 





x-1/2 > iV^ - 1. 



It follows from (A.2) and (A. 10) that 

< Rank(W*) < iVj - Rank(Z;;) < 1, 



(A. 10) 



(A. 11) 



that is, W* must be of rank one. 

What remains is to show that Y*„ is indeed of rank one. Firstly, one can show that c„ > since, if 
not, by (A.5) and by the fact of Q„„ >- 0, we must have Y*„ = 0, which, however, will lead to Z* ^ 
and W* = 0, and thus contradicts with (A.2). By substituting (A.8) into (A.4), we obtain the following 
two equalities 



— W: + A*„Q„„ ) A„ + —W^^KnK = 0, 

7n / 7n 

" n ~r '^nn^nn I On H »» n "■nnCn — " • 

In J In 



(A. 12) 
(A. 13) 



Further right multiplying (A. 13) with — b^/c„, and adding the resultant equality to (A. 12) gives rise to 

T— W; + K^Clnr)j ( A„ - hnh'^/Cn) = 0. (A. 14) 

Since (::^W* + K,n^nn) >- due to both A*„ > and Qn„ >- 0, (A. 14) implies that A„ = b„b^/c„. 



and thus 
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(A. 15) 
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which is a rank-one matrix. Case CI) is thus proved. 

Proof of case C2): Case C2) can be proved following similar derivations in (A. 10) and (A. 11), but 
using the KKT conditions of problem (16) with Qnnk = ooIaTj (i.e., e.„„fc = 0) for all n, k. 

Proof of case C3): Case C3) is a generalization of the result in [29] where the tightness of SDR for the 
worst-case robust beamforming problem in the single-cell scenario (Nc = 1) is studied. The condition in 
(17) can be proved following exactly the same idea as in [29] and thus the details are omitted here. ■ 
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